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Abstract 

A  method  to  calculate  the  three  principal  growth  stresses  in  Si02  scales  formed  during  SiC  fiber 
oxidation  has  been  developed.  The  method  assumes  that  during  oxidation  the  initial  volume  expansion 
at  the  SiC-Si02  interface  is  three-dimensional  and  equal  in  all  directions,  and  that  subsequent  Si02  shear 
stress  relaxation  is  described  by  the  stress-dependent  Eyring  viscosity  model.  Large  compressive  stresses 
of  ~10  GPa  in  Si02  adjacent  to  the  SiC-Si02  interface  are  relaxed  to  much  lower  levels  at  all 
temperatures.  At  1200°  -  1300°C  viscous  flow  of  amorphous  Si02  further  relaxes  stress  to  negligible 
levels.  At  700°  -  900°C,  axial  and  hoop  stress  at  the  GPa  level  persist.  Radial  stresses  only  reach  values 
greater  than  100  MPa  at  700°  -  900°C  for  scales  thicker  than  ~0.1  fiber  radii.  Radial  expansion  of  the 
scale  eventually  causes  hoop  stress  and  later  axial  stress  to  become  tensile  in  the  outer  scale. 
Differences  in  stress-states  developed  for  crystallized  and  uncrystallized  scales  are  considered.  Some 
tentative  calculations  for  crystalline  Si02  scales  are  compared  with  experimental  evidence  for  stress  in 
the  crystalline  Si02  scales  of  Hi-Nicalon™-S  SiC  fibers.  Assumptions  and  limitations  of  the  method  are 
discussed,  along  with  implications  for  fiber  strength  and  oxidation  kinetics. 

I.  Introduction 

Passive  oxidation  of  SiC  to  Si02  is  in  most  cases  controlled  by  interstitial  diffusion  of  02 
molecules  through  the  amorphous  or  crystalline  silica  scale.1"2  There  is  a  volume  expansion  of  ~2.2x 
during  SiC  oxidation  to  Si02.  Constraint  of  this  expansion  at  the  SiC-Si02  interface  generates  very  large 
stresses.  Microstructural  evidence  for  these  stresses  exists  for  crystalline  scales.  Cracks  formed  in 
crystalline  scales  during  SiC  fiber  oxidation  were  attributed  to  tensile  hoop  growth  stress  in  the  outer 
scale.3  5  These  cracks  function  as  short-circuit  diffusion  pathways  for  oxidation  and  render  the  scale  non¬ 
passivating.5  6  Tensile  stress  in  the  scale  may  also  contribute  to  lower  strengths  during  oxidation  at  high 
temperature.7  Very  high  dislocation  densities  in  crystalline  Si02  near  the  SiC-Si02  interface  scales 
suggests  that  very  high  shear  stresses  exist  during  scale  growth  near  that  interface.5 

Oxidation  of  silicon  to  Si02  has  a  volume  expansion  similar  to  that  for  SiC  oxidation.  Large  effort 
has  been  devoted  to  modeling  growth  stress  during  silicon  oxidation  at  700°  -  1200°C.8  These  studies  are 
motivated  by  the  use  of  Si02  as  a  dielectric  and  a  device  isolation  material  in  integrated  circuits.9'10 
Numerical  models  using  finite  element  or  finite  difference  methods  that  are  adaptable  to  a  variety  of 
silicon  substrate  geometries  in  integrated  circuits  have  been  developed  to  model  oxidation  growth 
stress.8'9, 11-19  Analytical  models  for  flat-plate  and  cylindrical  silicon  substrates  have  also  been  developed 
that  predict  differences  in  the  growth  stresses  for  the  two  geometries.  Radial  and  axial  compressive 
growth  stress  and  tensile  hoop  growth  stress  are  predicted  for  oxidation  of  silicon  fibers.20"22  For  the  flat 
plate  geometry,  only  in-plane  compressive  growth  stress  is  present. 

The  role  of  growth  stress  on  oxidation  mechanisms  in  silicon  has  been  extensively  analyzed  and 
discussed.23  Compressive  growth  stress  reduces  02  diffusivity  in  Si02,8,  12‘ 20,  24  25  reducing  Si  oxidation 
rates.26"28  Externally  applied  compressive  stress  has  also  been  shown  to  decrease  the  rate  of  silicon 
oxidation.29  Externally  applied  tensile  stress  increases  those  oxidation  rates,29'31  and  recently  the  same 
has  been  demonstrated  for  SiC  fibers.32  However,  growth  stress  modeling  and  analysis  for  SiC  oxidation 
is  lacking.  Extensive  characterization  and  analysis  of  carbon  and  SiOC  layer  formation  at  the  SiC-Si02 
interface  during  SiC  oxidation  has  not  considered  the  effects  of  growth  stress.2  Growth  stress  also 
contributes  to  the  residual  stress  in  Si02  scales  on  SiC  fibers,  which  in  turn  may  affect  the  crystallization 
rates  of  the  scales  and  the  strengths  of  oxidized  fibers.5 

Existing  growth  stress  models  for  cylindrical  silicon  substrates  assume  that  the  volume 
expansion  accompanying  oxidation  is  either  dilational  (equal  expansion  in  all  directions),7  uniaxial  (all 
expansion  in  the  radial  direction),7, 33-34  or  nearly  uniaxial  but  with  a  small  intrinsic  axial  and  hoop  strain 
in  the  Si02  oxidation  product.20, 25  Most  models  assume  stress  is  two-dimensional  and  ignore  the  axial 
stress  (ctz).21,  25, 33  34  However,  silicon  wafer  bending  during  oxidation  can  only  be  explained  by  three- 
dimensional  expansion  during  oxidation.8, 10  Early  cylinder  oxidation  models  assumed  Newtonian  viscous 
flow.7, 21, 33  Most  later  models  recognize  that  flow  at  high  stress  is  non-Newtonian  and  use  the  Eyring 
model  for  stress-dependent  viscosity.20,  25, 34  Uniaxial  expansion  models  with  Newtonian  flow  predict 
high  tensile  hoop  stress  (cre)  at  the  SiC-Si02  interface  that  decreases  towards  the  Si02  surface.7,  21,  33 
Models  with  stress-dependent  viscosity  predict  compressive  C7e  at  the  SiC-Si02  interface  that  rapidly 
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changes  to  tension  towards  the  Si02  surface.20, 34  One  model  with  stress-dependent  viscosity  considers 
only  shear  stress  and  does  not  attempt  to  resolve  it  into  principal  stresses.25  No  models  calculate  both  ctz 
and  crefor  cylindrical  substrates,  or  consider  how  the  change  in  c^from  compressive  to  tensile  affects 
shear  stress,  and  how  the  relaxation  of  shear  stress  then  separately  affects  <7eand  az .  For  structural  SiC 
fibers  this  is  an  important  omission  because  the  az  component  of  residual  stress  that  may  affect  fiber 
strength. 

A  method  for  computation  of  all  three  principal  growth  stresses  throughout  the  thickness  of  an 
Si02  scale  during  oxidation  of  SiC  fibers  is  presented.  The  method  can  be  equally  well  applied  to  silicon 
fibers  or  cylinders.  The  model  assumes  oxidation  volume  expansion  is  three-dimensional  (dilatational) 
and  that  the  stresses  resulting  from  constraint  of  that  expansion  are  relaxed  with  an  Eyring  stress- 
dependent  Si02  viscosity.35  Growth  stress  calculations  for  oxidation  of  SiC  fibers  are  done  from  700°  to 
1300°C  for  amorphous  scale  thicknesses  up  to  3  pm.  Some  calculations  are  also  done  for  crystalline 
scales  using  the  very  limited  data  available  for  cristobalite  viscosity.36"38  These  calculations  are  compared 
with  microstructural  evidence  for  stress  in  crystalline  scales  on  Hi-Nicalon™-S  SiC  fibers  that  has  been 
previously  published.5  Possible  effects  of  self  stress  on  02  diffusion  are  discussed.  The  change  in 
oxidation  growth  stress  with  fiber  radii  from  3  pm  to  ->  (flat  plate)  are  examined.  Assumptions  and 
limitations  of  the  model  are  discussed. 

II.  Method 

A.  General 

A  schematic  illustrating  the  volume  expansion  and  the  stresses  associated  with  oxidation,  as 
well  as  the  discretization  of  SiC  oxidation  into  annular  Si02  layers  is  shown  in  figure  1.  There  are  two 
sources  of  stress.  The  first  is  the  2.2x  volume  expansion  from  oxidation  of  SiC  to  Si02.  This  expansion  is 
shown  in  stress-free  and  constrained  states  in  figure  1.  Constraint  of  that  expansion  by  the  underlying 
SiC  and  overlying  Si02  causes  extremely  large  axial  (crz)  and  hoop  (ae)  stress,  and  much  smaller  radial  (<?r) 
stress  because  <7r  =  0  at  the  scale  surface  and  small  scale  thicknesses  relative  to  fiber  radii  are  usually  of 
interest.  crr  increases  or  decreases  through  the  scale  thickness,  depending  on  the  sign  and  magnitude  of 
ctz  and  ere.  The  stresses  involved  in  constraint  of  the  initial  volume  expansion  are  far  larger  than  those  for 
which  linear  elasticity  theory  is  valid.8  However,  as  pointed  out  previously,25  even  at  low  temperatures 
<7Z  and  ere  relax  very  rapidly  from  high  levels  by  radial  expansion  to  values  appropriate  for  linear  elasticity 
if  the  Eyring  model  for  the  dependence  of  viscosity  on  shear  stress  is  used.  The  second  source  of  stress  is 
the  circumferential  expansion  of  old  scale  as  it  is  radially  displaced  by  the  formation  of  new  scale  (Fig.  1). 
This  creates  tensile  <7@  in  the  older  scale  and  eventually  tensile  crz  by  the  Poisson  effect. 

A  flowchart  with  the  major  calculation  steps  is  shown  in  figure  2.  A  detailed  description  of  the 
calculation  methods  in  each  step  follows. 

B.  SiC  Fiber  Oxidation  Kinetics 

Fiber  oxidation  kinetics  do  not  deviate  significantly  from  those  for  flat  plate  geometry  until  the 
oxidation  products  for  12  pm  diameter  fibers  are  several  microns  thick,5,39  so  Si02  thickness  [w(t)j  (Fig. 
1)  is  assumed  to  obey  Deal-Grove  kinetics  for  flat  plate  geometry:5, 40 


w(t)  =  -A  1  + 


4B(t  +  ti) 


A2 


ti  = 


xf  +  Ax; 

B 


[la] 

[lb] 


where  ti  is  a  time  shift  that  corrects  for  the  presence  of  an  initial  oxide  layer  Xj  ,  B  is  the  parabolic  rate 
constant  and  B/A  is  the  linear  rate  constant.  All  growth  stress  modeling  for  amorphous  Si02  assumes 
that  X|=0,  but  for  crystalline  scales  it  is  useful  to  consider  the  thickness  at  which  the  scale  crystallizes  as  a 
stress-free  state  with  thickness  x,.  B  and  A  obey  the  usual  Arrhenius  relationships: 


-Qa 

A  =  A0e  RT 


[2] 


2 

Approved  for  public  release;  distribution  unlimited. 


[3] 


-Qb 

B  =  B0e  rt 


Fig.  1  Schematic  diagrams  of  oxidation  of  SiC  to  Si02.  Discretization  of  annular  layers  is  illustrated  and 
notated.  Volume  expansions  and  stresses  associated  with  oxidation  are  indicated.  The  circumferential 
expansion  of  an  old  scale  layer  as  it  is  pushed  outward  is  also  shown.  The  original  SiC  surface  is  shown  as 
a  blue  line,  roughly  midway  through  the  scale  thickness.  Original  SiC  is  dark  green,  the  elastically 
constrained  Si02  oxidation  product  is  light  red,  and  the  relaxed  Si02  is  light  green  with  a  purple 
boundary. 
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where  T  is  absolute  temperature,  R  is  the  gas 
constant,  Qa  and  QB  are  activation  energies  and  A0 
and  B0  are  pre-exponential  factors.  Recent  work  finds 
that  A0  =  4.7  x  10'4  m,  QA  =  110  kJ/mol,  Bp=  1.2  x  10'8 
m2/s,  Qb  =  248  kJ/mol  for  Hi-NicalonT  -S  SiC  fiber 
oxidation.5  The  Si02  oxidation  product  has  2.2  x  the 
volume  of  the  SiC  consumed.  The  SiC  radius  [b(t)] 
after  oxidation  for  time  t  is  (Fig.  1): 


b(t)  =  w2(n2  —  H)  +  —  flw 


[4] 


where  b0  is  the  original  fiber  radius  (6.1  pm  for  Hi- 
Nicalon™-S)5  and  O  is  the  ratio  of  SiC  molar  volume 
(QSic)  to  Si02  molar  volume  (fiSi02b  sometimes  referred 
to  as  the  Pilling-Bedworth  ratio. 


n  = 


flc; 


SiC 


12.46  cm3 


=  0.456 


^sio2  27.34  cm3 
The  outer  radius  of  the  Si02  scale  [c(t)]  is  (Fig.  1): 

c(t)=b(t)+w(t) 


[5] 


[6] 


C.  Elastic  Growth  Stress 

The  elastic  stresses  and  strains  can  be 
determined  from  a  modification  of  a  method  used  for 
progressively  deposited  coatings  on  cylinders.41  This 
method  was  developed  for  sequentially  deposited 
coatings,  where  the  last  coating  is  at  the  surface. 
During  fiber  oxidation  the  "deposition"  sequence  is 
reversed;  the  last  oxide  increment  forms  at  the  SiC- 
Si02  interface.  The  Si02  scale  is  therefore  discretized 
into  annular  increments,  or  layers,  formed  in  equal 
units  of  time,  with  the  oldest  layer  at  the  surface. 

The  unrelaxed  elastic  growth  stresses  are 
calculated  by  dividing  the  SiC-Si02  system  into  three 
separate  volumes  -  the  unoxidized  SiC  fiber,  the 
annular  Si02  added  from  time  t(i  -1)  to  t(i)  at  the  SiC- 
Si02  interface,  and  the  outer  (older)  Si02  ring  that  was 
added  from  time  t(0)  to  t(i-l)  (Fig.  1).  The  axial,  radial, 
and  hoop  stresses  in  the  SiC  fiber  at  time  t(i)  are 
designated  crzSlC(i),  arSlC(i)  .,  and  qnS'c(i)  respectively. 
The  axial,  radial,  and  hoop  stresses  in  the  two  annular  Si02  layers  are  crzSl02(i),  <jrSl  2(i)  ,  ae5l02(i)  and 
CTzSl02(i-l),  crrSl02(i-l) ,  cJeSl02(i-l),  respectively.  The  increment  in  Si02  thickness  formed  from  time  t=i-l  to 
t=i  is  therefore  w(i)-w(i-l).  The  effect  of  adding  the  ith  Si02  layer  on  the  principal  stresses  in  the  three 
separate  volumes  can  be  found  from  three  equations  for  f,  Pj,  and  pis,  using  a  modification  of  the 
method  of  Tsui  and  Clyne  that  uses  strain  compatibility  equations:41 


growth  stresses  during  SiC  fiber  oxidation. 


ESio2  =  a2‘°2  (i  -  1)  -  g2‘°2  (i)  +  uSi02  (ag‘°2  (i)  -  Gg‘°2  (i  -  1)  +  af°2  (i)  -  a^2  (i  —  1))  [7] 

ESio2 Ae  =  Oq °2 0  -  1)  -  oSq°2 (0  +  uSio2  (^z‘°2 (0  -  °z°2 0  -  1)  +  r‘°2 (0  -  r‘°2 0  -  1))  [8] 
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—  [^e°2(i)  -^Si02(^zl0z(i)  +  tfrlO20))]  =  [®|lC(0  -  USic(<^|1C  (0  +  ^riC(i))] 


where: 


.  3  ^Si02 

As  =  — — --  1 

^ SiC 


°ziC  0)  - 


-Esicf 


u  ^ESicb2(i)  +  Esio2  ((b(i)  +  w(i)  -  w(i  -  l))2  -  b2(i)) 

CTriC0)  =  crfiC(i  —  1)  —  PiS 

GelC0)  =  a|lC(i-l)-pis 

Si02r.x  _ESi02f 

°z  W=-7 - 2 - \ 

n  ^ESicb(i)2  +  ESiQ2  ((b(i)  +  w(i))  -  b(i)2)  j 


+  <jfic(i-  1) 


Si02  /■■  /<  \ 

Gz  2(l  —  1)  = 


Si02 

Gr  (0  =  “Pi 


siozn  _  Pj(b(i)  +  w(0) 
Ge  uj  w(i)  -  w(i  -  1) 

f 


Si02  /'•  y* 

Gfi  2  (l  -  1)  = 


tt  (b(i  —  1)  +  w(i  —  1)  —  w(i))2  —  b2(i  —  1)| 


Gri02(i-1)  =  0 


2pisb(i  -  l)2  -  Pi  ((b (i  -  1)  +  w(i  -  l))2  +  b(i  -  l)2) 
(b(i  -  1)  +  w(i  -  l))2  -b(i-l)2 


[9] 

[10] 

[11] 

[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 


where  f  is  the  axial  force,  pis  is  the  pressure  across  the  Si02-SiC  interface,  p,  is  the  pressure  across  the 
interface  between  the  ith  and  the  i-lth  Si02  layers,  ESiC  and  ESi02  and  Young's  modulus  of  the  SiC  fiber  and 
the  Si02  scale,  respectively,  and  vSic  and  vSi02  are  Poison's  ratio  for  the  SiC  fiber  and  the  Si02  scale. 
Stresses  in  older  layers  (j  =  i-2  to  j  =  0)  are  updated  with  the  i-1  stress  values  in  [17  - 19]: 

Gzi02(j)  =  a*i0’(j)+a?°’(i-l)  [20] 

af°2(j)  =  af°2a)+  0riO2(i-l)  [21] 

aSei02(i)  =  G0°2(j)  +  Gei02(i-D  [22] 


D.  Relaxation  of  Shear  Stress 

The  relaxation  of  the  elastic  stresses  for  all  annular  elements  (j=l  to  i)  in  time  increment  At=t(i)- 
t(i-l)  are  calculated  next.  Stress  relaxation  in  the  SiC  substrate  is  assumed  to  be  negligible.  Deviation 
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from  Newtonian  viscosity,  referred  to  as  glass  plasticity,  results  in  lower  effective  viscosity  at  high  shear 
stress.  The  Eyring  model  for  shear  stress  (t)  dependence  of  glass  viscosity  (ri)  is  frequently  used  for 

viscosity8"9' 11-*2, 14, 16, 18-20, 22, 25, 35, 42 


8  =  Tlo 


Sinh(T/xc) 


[23] 


where  the  activation  volume  for  plasticity  in  Si02,  decreases  with  temperature,  and  has  been 
inferred  to  have  values  ranging  between  1.2  xlCP8  rn  to  3  x  10'28  m3.12,  20,  25  The  parameter  k  is 
Boltzmann's  constant,  xc  is  the  critical  shear  stress  above  which  plasticity  is  significant  (typically  ~  150 
MPa),  and  r|0  is  the  stress-fee  viscosity  described  by:43"44 

Q 

r]0  =  C0  eRT  ,  C0  =  3.8  x  1(T13  &  Q  =  712,000  kj/mol  [24] 

from  1000°C  to  1400°C.  Unfortunately  this  relationship  lacks  experimental  validation  below  1000°C. 43-44 
For  silica  formed  from  SiC  oxidation,  carbon  in  the  Si02  scale  may  cause  viscosities  to  be  higher  than 
those  in  [24], 45  The  plasticity  activation  volume  Vc  corresponds  to  a  critical  shear  stress  xcof  about  100 
MPa,  which  is  roughly  consistent  with  experiment  at  500  to  1400°C.42'  46  Shear  stress  relaxation  with 
time  obeys  a  Maxwell  viscoelastic  model:12, 25,47 


dx(t)  _  — G  x(t)  / 
dt  AlCO  ' 


x[t(i  -  1)]  =  x0  t(i  -  1) 


[25] 


where  G  is  the  Si02  shear  modulus  and  the  shear  stress  at  the  initial  time  t(i-l)  is  t„.  The  shear  modulus 
(G)  of  amorphous  silica  is  34  GPa,  with  only  weak  temperature  dependence.48  The  relaxation  of  shear 
stress  with  time  (x(t))  can  be  determined  by  substitution  of  [23-24]  in  [25]  and  solving  the  resulting 
differential  equation  [25]  for  x[t]: 


The  initial  shear  stress  x0  is  determined  from  the  principal  stresses  for  all  the  separate  annular  elements 
(j=l  to  i)  by  the  usual  method: 

x0  =^Oi(j)  -  o3(j))  [27] 

where  cri(j),  a2(j),  and  a3(j)  are  the  maximum,  intermediate,  and  minimum  principal  stresses  of  the  jth 
annular  element  at  time  t(j)  determined  from  values  for  azSl02(j),  CTrSl02(j),  and  aeSl02(j). 

E.  Determination  of  Relaxed  Principal  Stresses  from  Relaxed  Shear  Stress 

Calculation  of  the  relaxation  of  the  principal  stresses  from  relaxation  of  the  shear  stress  x0  to  x(t) 
in  equation  [26]  is  complicated  by  the  presence  of  both  tensile  and  compressive  principal  stresses.  It  is 
therefore  also  useful  to  designate  principal  stresses  by  their  rank  in  absolute  value  as  (j) ,  2a  (j),  and 
3a(j),  and  to  define: 


s[o®i={-i^<S  1281 

A  boundary  condition  is  arSl02(l)  =  0,  but  arSl02(j>l)  will  have  a  small  negative  value 
(compression)  if  cjzSl  2(j>l)  and  aeSl02(j>l)  are  positive  (tension),  or  a  positive  value  if  the  other  principal 
stresses  are  negative.  Since  the  scale  is  free  to  expand  in  the  radial  direction,  volume  is  not  conserved, 
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and  decomposition  of  the  stress  tensor  into  hydrostatic  and  deviatoric  components,  as  used  for  some 
growth  stress  calculations,13'14  is  not  useful  because  the  hydrostatic  pressure  can  change  with  relaxation. 
In  almost  all  cases  arSl02(j)  will  be  3a  (j),  but  there  are  some  exceptions  when  gzs'°  (j)  or  aeSl02(j)  are 
changing  from  compressive  (negative)  to  tensile  (positive)  and  are  therefore  3a  (j)  for  a  relatively  short 
period  of  time.  The  different  relationships  between  Gi,2,3(j)  and  1,2'3<j(j)  are  shown  in  figure  3,  and 
require  differences  in  the  way  new  principal  stresses  [a(j)'j  are  calculated  after  shear  stress  relaxation. 
Two  methods  are  employed.  The  first  assumes  that  only  3aQ)  will  relax,  until  it  reaches  the  value  for 
2a(j),  at  which  point  both  relax.  This  method  is  best  suited  for  high  x  where  q  [23]  is  not  Newtonian,  and 
is  later  referenced  as  the  "High  x  Method".  arSl02(j)  is  a  boundary  condition  unless  |arSl02(j)|  <  |azSl02(j)| 
&  |CTfiSl02(j)  |,  where  the  boundary  condition  becomes  0  for  all  stresses.  For  principal  stresses  of  type  (Sji, 
2a2,  3g3)  or  (3CTi,  2g2,  3g3): 


If  2x(j)  >  |2aO)-3a(j)  | : 

If  2x(j)  <  1 2cr(j)  -  3a(j)  |  &  8[2cr(j)]  =  S[3a(j)]: 

If  2xQ)  <  I2a(j)  -  aG)l  &  S[2a(j)]  *  S[3gG)]  &  x(j)  > 
If  2xQ)  <  1 2a(j)  -  3^G)  I  &  S[2a())]  *  8[3g(])]  &  x(j)  < 


For  all  other  principal  stresses: 

lfx(j)>  I  2ctQ)  I : 
if  l2c>G)l  >x(j)>  l3c>G)|: 
If  tQ)  <  I  3ctG)  I : 


ctG)'  =>G)  +2S[1gG)]x 

°0)'=  °G) 

3a())'  =  3a0 

[29] 

W  =  >G)  +  28[1aG)]i 

W  =  >G)  +  2S[1a(j)]x 

3a(j)'  =  3aG) 

[30] 

>G)'  =  >G)  +  28[1a(j)]T 
2M})'=  ctG)  +  28[1gG)]x 

3aG)'  =  ctG) 

[31] 

>G)'  =  5[Wj)]T 

aG)'  =  8[  G(j)]x 

3G(j)'  =  S[3G(j)]T 

[32] 

3G(j)'  =  2gG)  +  28[1G(j)]t 
aG)'=  aG) 

3a(j)'  =  3aG) 

[33] 

1G(j)'  =  8[1G(j)]x 
gG)'  =  8[2gG)]x 

3aG)'  =  3aG) 

[34] 

>G)'  =  8[!g(])]x 
a G)'  =  S[  G(j)]x 
aG)'  =  S[3g(])]x 

[35] 

This  method  is  consistent  with  radial  expansion  at  a  free  surface  with  minimization  of  strain  energy 
through  minimization  of  3gG)2  +  2gQ)2  +  3gG)2-  The  relaxed  values  of  ^(j)',  2cr(j)',  and  3ctQ)'  are 
reassigned  to  azSl02(j),  arSl02(j),  and  aes'°2(j). 

The  second  method  assumes  that  viscosity  is  Newtonian  and  that  the  rate  of  Ge(j)and  az(j) 
relaxation  is  proportional  to  their  difference  with  x(t)/x0and  Gr(j)  ,  which  is  a  boundary  condition  and 
does  not  change,  being  zero  at  the  Si02  surface  and  near  zero  elsewhere.  This  method  is  best  suited  for 
high  temperatures  and  low  shear  stresses  where  q  is  Newtonian,  and  is  later  referenced  as  the  "Low  x 
Method",  ffe(j)'  and  az(j)'  can  be  determined  by  solution  of: 


a0(j)' 

0 

0 

_  T(t) 

a0(j)  -arQ) 

0 

0 

ar(j) 

0 

0 

0 

az(j)' 

0 

0 

CTz(j)  -  ar(j) 

0 

+ 

0 

CTr(j) 

0 

0 

0 

CTr(j). 

To 

0 

0 

GrG). 

0 

0 

ar(j). 
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At  very  low  principal  stresses,  where  |arSl02(j)|  >  JarzSl02(j)|  & 
I  CT0SlO2(j)  I  and  they  are  opposite  in  sign,  |az  (j)|  and 
I CT0SlO2(j)  |  increase,  and  this  method  is  therefore  inconsistent  with 
minimization  of  strain  energy.  However,  this  condition  can  only 
occur  very  briefly  when  azSl02(j)  and  creS'02(j)  change  from 
compressive  to  tensile  from  radial  displacement  of  the  scale, 
described  next. 

F.  Radial  Displacement  and  Hoop  Stress  Generation  in  Scale 

Relaxation  expands  the  Si02  scale  radially.  The  individual 
radial  displacement  of  the  jth  layer  (ur(j) )  is: 


ur(j)  = 

where 


n 


SiO-, 


nSiC  (l  +  ^i02G)  +  4i02G)  +  eli02(j)eei02(j)) 


-  1  [37] 


e|i02(j)  =  ^-kSi02(j)  -  uSio2  ki02(j)  +  arSi02(j))l  [38] 

bSi02 

4i02G)  =  F^ki02G)  -  ^Si02  kSi02G)  +  ^i02(j))l  [39] 

bSi02  L  v 


Fig.  3  The  relationships  between 
the  magnitudes  (crj)  and  absolute 
value  magnitudes  ('a)  of  the 
principal  growth  stresses  that  are 
used  to  calculate  principal  stress 
relaxation  from  shear  stress 
relaxation  in  [29]  -  [35], 


The  radial  displacement  of  the  jth  outer  layer  is  equal  to  the  sum  of 
the  radial  displacements  of  all  the  younger  layers  (i-j-1)  beneath  it. 
This  adds  hoop  strain  (seSl02)  to  outer,  older  layers  as  they  are 
displaced  radially  outward  and  forced  to  a  larger  circumference 
(Fig.  1).  The  new  hoop  strain  in  each  layer  can  be  calculated  from 
the  sum: 


where 


„Si02 


G)  = 


,Si02 


G)  + 


i 

^  ur(j) 


b(j  ~  1)  ~  b(j) 
b(i) 


w(i)  -  w(i  -  1)  =  ^  ur  (i)(b(i  -  1)  -  b(i)) 

i-l 


[40] 


[41] 


for  complete  relaxation  of  the  youngest  layer. 

G.  Recalculation  of  Elastic  Stress  in  the  Scale  and  SiC  Fiber  after  Radial  Displacement 

The  stress  in  each  Si02  layer  can  then  be  recalculated  from  the  new  eeSl02  by  solving  the  three 
strain  compatibility  equations  for  the  three  principal  stresses: 

tf'02(i)  =  F^-k'02®  -  uS,o,  k'02ffi  +  °?°2®)]  [42] 

hSi02  L  v  /J 

4i02  G)  =  F^-ki02G)  -  ^Si02  kSi02G)  +  ^i02G))l  [43] 

bSi02  L  v 

epi02(j)  =  F“ —  ki02G)  —  uSi02  (CT!iO20)  +  t*zi02®)l  [44] 

bSi02  L  v  /J 
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The  revised  principal  stresses  in  the  SiC  fiber  can  then  be  calculated.  The  axial  stress  azSlC  is  computed 
from  the  total  force  exerted  by  the  Si02  layers  which  is  the  sum  of  the  axial  stress  in  each  annular  Si02 
element  x  area  of  that  element: 


i-l 


CTziC0)  =  —  ^ 

j=i 


gzi02(j)[2c(i)(w(j)  -  w(j  -  1)  +  w(j)2  -  w(j  +  l)2] 

b(i)2 


[45] 


The  revised  radial  and  hoop  stress  in  the  SiC  fiber  can  be  computed  in  a  similar  manner,  by  computing 
the  net  pressure  (pn)  from  the  sum  of  the  pressures  in  each  annular  layer: 


CTriC(0  =  cr|lC(0  =  “Pn 


These  revised  stresses  are  changed  incrementally  for  strain  compatibility  and  stress  relaxation  as  the 
program  loops  back  to  equations  [1  -  41]  as  the  next  layer  i  is  considered.  For  SiC  oxidation  when  w(i) 
«  b(i),  the  stress  in  Si02  is  much  larger  than  that  in  SiC  and  can  be  neglected  for  approximate 
calculations. 


i-l 


—  I' 


_Si02 

re 


(i) 


j=i 


w(j)  -  w(j  -  1) 
b(i) 


[46] 


III.  Results  and  Discussion 
A.  General 

Growth  stress  calculations  were  done  using  equations  [1-46]  in  a  Mathematica™  program. 


Calculations  were  done  with  50  to  2000  layers  ( 
insignificant  for  i  >  100.  A  value  of  i  =  500  was 
calculation  speed. 

Calculations  were  done  for  oxidation  of  h 
1000,  and  3000  nm  thickness  (w)  at  700,  800, 


Distance  from  SiC-Si02  Interface  (nm) 

Fig.  4.  Floop  stress  (oe)  and  axial  stress  (ctz)  for 
a  1  pm  thick  scale  on  a  Hi-Nicalon™-S  (6  pm 
radius)  fiber  calculated  using  the  "High  x 
Method"  [29-35]  and  "Low  x  Method"  [36], 
Differences  in  calculated  values  between  the 
two  methods  are  small.  Compressive  values 
are  negative. 


i  =  50  to  2000).  Variations  in  calculations  were  usually 
chosen  as  an  optimal  trade-off  between  precision  and 

li-Nicalon™-S  SiC  fiber  for  amorphous  scales  of  10,  100, 
900,  1000,  1100,  1200,  and  1300°C.  The  Deal-Grove 
oxidation  kinetics  for  this  fiber  have  been  reported 
for  dry  air  between  700  and  1300°C  with  A0  =  6.5  x 
10'4  m,  B0  =  1.2  x  10'8  m2/s,  QA  =  111  kJ/mol,  and  QB  = 
249  kJ/mol.5  6  The  fiber  radius  (b0)  is  6.1  pm.  The 
molar  volumes  for  SiC  (nSiC)  and  amorphous  Si02 
(Clsim)  are  27.34  cm3  and  12.46  cm3,  respectively. 
The  Young's  modulus  (E)  and  Poisson's  ratio  (v) 
values  used  for  SiC  and  Si02  were  ESiC  =  400  GPa,  vSic  = 
0.157,  ESio2  =  70  GPa,  and  vSi02  =  0-17.  The  shear 
modulus  value  used  for  silica  (G)  was  34  GPa.  The 
continuous  variation  in  growth  stress  with  change  in 
fiber  radius  was  examined  by  calculations  using  b0  =  3 
pm  and  b0  =  4°°  (1  km). 

General  features  for  hoop  stress  (cie)  and  axial  stress 
(ctz)  are  evident  in  figure  4  for  a  1  pm  thick  scale  on  a 
Hi-Nicalon™-S  (6  pm  radius)  fiber.  Qualitatively, 
these  results  are  most  similar  to  those  of  Delph  for 
elastic-viscoplastic  scales  during  silicon  oxidation.20 
Compressive  elastic  stress  of  ~-25  GPa  for  ere  and  az  is 
rapidly  relaxed  by  the  Eyring  model  [23]  for  Si02  glass 

viscosity  (r|)  to  ~-3  GPa  at  700°C  and  -1  GPa  at 
1000°C  just  10  nm  away  from  the  SiC-Si02  interface. 
This  quickly  reduces  the  problem  to  one  similar  to 
unidirectional  radial  expansion  with  a  residual 
"intrinsic  strain"  parallel  to  the  SiC-Si02  interface. 
However,  this  "intrinsic  strain"  is  T  and  t  dependent, 
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and  the  calculations  done  using  unidirectional  radial  expansion  show  only  tensile  cr0,  with  no  transition 
from  compressive  to  tensile  stress. 20,25 

At  the  scale  surface,  cr0is  ~2  GPa  (tension)  at  700°  and  800°C,  and  crz  is  ~1  GPa  (tension)  at 
700°C.  ct2  is  driven  towards  tensile  values  by  the  Poisson  effect  from  radial  expansion  that  creates 

tensile  a0.  Viscous  stress  relaxation  reduces  all  growth 
stresses  to  <  100  MPa  at  1100°C  and  <  10  MPa  at 
1300°C  everywhere  in  the  Si02  scale  except  the  region 
immediately  adjacent  to  the  SiC-Si02  interface. 

Radial  stress  (ar)  calculated  by  the  "High  x 
Method"  [29-35]  for  a  1  pm  thick  scale  on  a  Hi- 
Nicalon™-S  (6  pm  radius)  fiber  is  shown  in  figure  5 
and  shear  stress  (x)is  shown  in  figure  6.  Compressive 
CTr  up  to  ~100  MPa  is  present  in  the  center  of  the  scale 
at  700°  -  900°C.  This  is  a  consequence  of  tensile  cr0  at 
the  scale  surface.  At  700°  -  800°C  crr  goes  tensile  at 
the  SiC-Si02  interface,  which  is  a  consequence  of  the 
high  compressive  cr0  near  that  interface. 

x  of  ~12  GPa  is  rapidly  relaxed  to  ~  2  GPa  at 
700°C  and  500  MPa  at  1100°C  just  10  nm  away  from 
the  SiC-Si02  interface  (Fig.  6).  It  continues  to  decrease 
as  compressive  cr2  and  a0 decreases  away  from  the 
SiC-Si02  interface  until  cj0  becomes  tensile,  at  which 
point  x  increases  towards  the  Si02  surface,  reaching 
surface  values  close  to  1.5  GPa  at  700°C  and  300  MPa 
at  1000°C.  The  shear  stress  minimum  is  ~  600  nm 
from  the  SiC-Si02  interface  at  700°C,  200  nm  at 
1000°C,  and  ~50  nm  at  1100°C.  Note  that  a  1  pm 


Distance  from  SiC-Si02  Interface  (nm) 


Fig.  5.  Radial  stress  (ar)  for  a  1  pm  thick  scale 
on  a  Hi-Nicalon™-S  (6  pm  radius)  fiber 
calculated  using  the  "High  x  Method"  [29-35], 


Distance  from  SiC-Si02  Interface  (nm) 


scale  takes  ~109  s  to  form  at  700°C  and  ~106  s  to  form 
at  1000°C,  so  the  latter  is  of  more  practical  interest. 

The  practical  use  of  the  growth  stress 
calculation  method  is  likely  to  be  limited  by 
knowledge  of  accurate  values  for  amorphous  silica 
viscosity.  Unlike  silicon  oxidation,  silica  scales  formed 
during  SiC  oxidation  incorporate  carbon.49'52  Network 
carbon  in  amorphous  Si02  stiffens  the  network 
structure,  making  it  more  viscous  and  less  permeable 
to  02.45, 53  54  However,  Si02  viscosity  may  be  reduced 
by  incorporation  of  the  Cl  (620  ppm),  S  (52  ppm),  Ca 
(45  ppm),  Na  (35  ppm),  and  Fe  (27  ppm)  impurities  in 
Hi-Nicalon™-S  fiber.6 

B.  Comparison  of  Methods  for  Calculation  of 
Relaxed  Principal  Stresses 

Two  methods  for  calculation  of  the  relaxed 
principal  stresses  from  the  relaxed  shear  stress  were 
developed  (equations  [29-36]).  These  methods  are 
compared  for  a  1000  nm  thick  amorphous  Si02  scale 
on  a  6  pm  radius  Hi-Nicalon™-S  fiber  at  700°  -  1300°C 
in  figure  4.  The  differences  between  the  methods  are 
illustrated  by  the  thickness  of  the  colored  lines.  For  a0 
there  is  little  difference  between  the  two  methods. 
The  "Low  x  Method"  [36]  finds  cr0  values  slightly  lower 
[29-35],  In  the  compressive  stress  regions  of  the  scale 


Fig.  6.  Shear  stress  (x)  for  a  1  pm  thick  scale 
on  a  Hi-Nicalon™-S  (6  pm  radius)  fiber 
calculated  using  the  "High  x  Method"  [29-35], 

in  absolute  value  than  the  "High  x  Method 
(negative  values,  close  to  the  SiC-Si02  interface)  the  calculated  values  for  axial  stress  (az)  are  nearly 
identical.  Close  to  the  surface  it  is  higher  (more  tensile)  at  700°  and  800°C  and  lower  (less  tensile)  for 
higher  temperatures.  At  1100°C  and  higher,  stresses  are  rapidly  relaxed  by  viscous  flow.  There  is  a  2x 
relative  difference  in  calculated  values  of  <rz  between  the  two  methods,  but  the  absolute  difference  is 
small  because  stresses  are  low  (~  100  MPa)  and  highly  relaxed  by  viscous  flow.  Since  high  scale  stresses 
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are  generally  of  more  interest  for  structural  SiC  fibers,  further  calculations  are  presented  from  the  “High 
x  Method". 


C.  Variation  with  Scale  Thickness  and  Fiber  Radius 

Calculations  for  3  pm  radius  (b0),  6  pm  radius,  and  a  flat  plate  for  10  nm,  100  nm,  and  1  pm  thick 
(w)  scales  at  700  -  1300°C  are  shown  in  figure  7.  The  continuous  change  of  a0and  crz  with  change  in  b0 


Distance  from  SiC-Si02  Interface  (nm)  Distance  from  SiC-Si02  Interface  (nm)  Distance  from  SiC-Si02  Interface  (nm) 

Fig.  7.  Hoop  (ae)  and  axial  oxidation  growth  stress  (ctz)  levels  as  a  function  of  position  in  the  Si02  scale 
for  SiC  flat  plates  and  fibers  with  6  and  3  pm  radius  (b0).  Results  for  w  =  10  nm,  100  nm,  and  1  pm 
thick  scales  are  shown  for  temperatures  of  700°,  800°,  900°,  1000°,  1100°,  1200°,  and  1300°C.  The 
continuous  change  in  both  ere  and  ctz  with  the  variables  of  scale  position,  temperature  (T),  scale 
thickness  (w),  and  fiber  radius  (b0)  is  evident. 

and  w  at  700°C  to  1300°C  throughout  the  scale  is  evident.  For  flat  plates  of  CTe=  ctz,  and  stress  is  always 
compressive.  Growth  stress  for  w  =  10  nm  is  very  similar  to  that  of  a  flat  plate  for  b0  =  3  and  6  pm. 
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For  a  3  pm  radius  fiber,  tensile  a0  develops  at  the  scale  surface  at  temperatures  over  950°C  for 
w  =  100  nm  and  reaches  values  of  400  MPa  at  1100°C.  Tensile  cyz  does  not  develop  until  T  >  1150°C  and 
reaches  100  MPa  at  1200°C.  For  a  6  pm  radius  fiber,  tensile  stress  does  not  develop  at  the  scale  surface 
until  temperatures  close  to  1200°C  are  reached;  maximum  tensile  cr0and  az  of  100  MPa  forms  at 
1200°C.  For  both  b0  =  3  and  6  pm  with  w  =  100  nm,  crz  is  generally  much  more  compressive  than  a0, 
except  at  the  highest  temperatures  stresses  are  almost  completely  relaxed. 

The  results  for  w  =  1  pm  scales  with  b0  =  6  pm  were  discussed  in  the  previous  section.  Maximum 
tensile  ae  reaches  2.3  GPa  and  3  GPa  for  b0  =  6  and  3  pm,  respectively,  at  700°C.  Maximum  tensile  cz 
reaches  1  GPa  (900°C)  and  2.2  GPa  (800°C)  for  b0  =  6  and  3  pm,  respectively.  Much  larger  fractions  of  the 
scale  are  under  tensile  stress  for  a  1  pm  thick  scale  in  comparison  to  a  100  nm  thick  scale. 


D.  Surface  Stress 

For  structural  SiC  fibers  the  residual  stress  at  the  scale  surface  is  of  great  interest.  Tensile  stress 
promotes  formation  of  cracks  that  can  degrade  fiber  strength.5'6  Compressive  stress  may  have  the 
opposite  effect.  The  stress  values  at  maximum  distance  from  the  SiC-Si02  interface  in  figure  7  are  the 
surface  stresses  (right  hand  side  of  each  plot).  Calculations  for  the  hoop  stress  (ce)  and  axial  stress  (ctz) 
at  the  scale  surface  are  mapped  as  a  function  of  temperature  (T)  and  scale  thickness  (w)  in  figures  8-10 
<7Z  -  6  pm  Radius  _  Gz  -  3  pm  Radius 
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Fig.  8.  Axial  stress  (ctz)  at  the  Si02  scale  surface  as  a  function  of  scale  thickness  and  oxidation 
temperature  for  6  and  3  pm  radius  SiC  fibers.  Oxidation  time  isochrones  are  also  shown. 


Fig.  9.  Hoop  stress  (cr0)  at  the  Si02  scale  surface  as  a  function  of  scale  thickness  and  oxidation 
temperature  for  6  and  3  pm  radius  SiC  fibers.  Oxidation  time  isochrones  are  also  shown. 


12 

Approved  for  public  release;  distribution  unlimited. 


for  flat  plates,  6  pm  radius  fibers,  and  3  pm  radius  fibers.  These  plots  clearly  distinguish  tensile  and 
compressive  regions,  and  also  show  logarithmic  time  isochrones  in  seconds. 

The  tensile  stress  maximum  forms  at  smaller  scale  thicknesses  as  temperature  increases.  This 
maximum  is  the  greatest  practical  concern  for  structural  SiC  fibers.  The  plots  illustrate  T-t  windows  that 
should  be  avoided.  For  example,  a  tensile  crz  of  ~500  MPa  forms  for  a  scale  of  ~  600  nm  thickness  at 
1050°C  in  105  seconds  (~  1  day)  during  oxidation  of  a  6  pm  radius  fiber  in  dry  air  (Fig.  8).  For  a  3  pm 
radius  fiber  this  maximum  is  ~750  MPa.  Longer  oxidation  times  would  decrease,  rather  than  increase, 
this  stress.  Much  higher  tensile  stresses  are  possible  at  lower  temperatures,  but  require  much  longer  to 
form,  and  will  form  for  thicker  rather  than  thinner  scales.  These  high  tensile  stresses  may  impact  fiber 
strength  at  temperature.  They  are  larger  than  the  compressive  thermal  stresses  that  will  develop  in  the 
Si02  scale  during  cool  down  (~300  MPa),  and  so  may  impact  room  temperature  strength  as  well. 

The  maximum  tensile  a0  that  forms  at  1050°C  is  ~600  and  900  MPa  for  6  and  3  pm  radius  fibers, 
respectively  (Fig.  9).  This  stress  maximum  occurs  at  a  significantly  lower  scale  thickness  and  shorter  time 

than  the  crz  maximum  (Fig.  8).  Tensile  (?e 


3000- 


"a„&az"-  Flat  Plate 


Fig.  10.  Stress  (  ae  =  crz)  at  the  Si02  scale  surface 
as  a  function  of  scale  thickness  and  oxidation 
temperature.  Oxidation  time  isochrones  are 
also  shown. 

oxidation  volume  increase  is: 


stresses  have  been  suggested  to  cause  axial 
cracking  of  scale  during  oxidation,  rendering  the 
scale  non-passivating.  6  Calculations  for  a  flat 
plate,  where  all  stresses  are  compressive,  are 
shown  in  figure  10  for  comparison. 

E.  Steady-State  Tensile  Stress 

A  near  "steady-state"  develops  for 
tensile  stress  at  T  >  1200°C  for  w  =  100  nm  for 
both  b0  =  3  and  6  pm  (Fig.  7).  At  lower 
temperatures,  greater  scale  thickness  is  required 
to  reach  steady-state.  A  w  >  ~1  pm  is  required 
for  T  >  1000°C.  The  tensile  cr9and  crz  values  are 
nearly  equal,  with  little  change  for  significant 
distances  beneath  the  scale  surface  (Fig.  7),  but 
the  scale  thickness  (time)  required  to  reach 
steady-state  is  longer  for  ctz  because  it  is  driven 
to  tensile  values  by  the  Poisson  effect  from  ct9. 
For  shear  stress  x  <  tc  of  100  MPa  [23]  where  a 
stress-free  viscosity  (r|0)  [24]  is  applicable,  a 
simple  expression  for  the  steady-state  hoop 
stress  [ae(ss)]  can  be  derived.  The  rate  at  which 
scale  thickness  (5)  is  increased  from  the 


dt  w  Vn  / 


[47] 


If  <yr  is  negligible,  the  shear  stress  is  half  the  hoop  stress,  and  the  rate  at  which  hoop  strain  (se)  develops 
in  the  outer  scale  is: 


dfe 

dt 


x 

*lo 


aa(ss) 

2rio 


d6 

dt 

b 


[48] 


[47]  and  [48]  can  be  solved  for  the  steady-state  hoop  stress  g0(ss): 

2B(n  ~  i)i« 


Ge  (ss)  = 


wb 


[49] 
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For  x  >  100  MPa  the  stress  dependence  of  viscosity  is  significant,  and  [23]  must  be  substituted  for  r|0, 
giving: 


Ge  (ss)  =  — 2xcCsch 


bwxca 

.Bti0(n-  1). 


[50] 
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Predictions  for  steady-state  ae tensile  stress  from 
[50]  for  a  6  pm  radius  fiber  are  shown  in  figure  11. 
High  temperature  (>  1150°C)  steady-state  stress 
value  predictions  are  high,  but  at  lower 
temperatures  the  prediction  is  close  to  calculated 
values.  A  steady  state  develops  in  ~2  minutes  at 
1300°C,  3  hours  at  1100°C,  and  in  ~100  days  at 
900°C.  The  predictions  for  Gzare  the  same  as  those 
for  ae,  except  the  times  to  reach  steady  state  are 
longer  (Fig.  8).  For  smaller  radius  fibers  the  time  to 
reach  steady-state  is  also  shorter. 

By  expanding  [49]  with  the  Arrhenius 
expressions  for  B  and  r|0,  the  steady-state  tensile 
stress  decrease  with  increasing  temperature  is 
evident  when  Q  >  Qb;  the  scale  viscosity  decreases 
faster  with  an  increase  in  temperature  than  the 
oxidation  rate: 


G0  (ss)  = 


Q-Qb 
2B0C0e  RT 


wb 


[51] 


Fig.  11.  The  plot  for  surface  stress  in  Fig.  9  for  b0 
=  6  pm,  with  predictions  for  the  steady-state 
region  from  [50]  superimposed. 

generation  associated  with  it,  will  increase  faster 


This  suggests  that  if  materials  exist  for  which  Q  < 
Qb,  a  counterintuitive  increase  in  tensile  growth 
stress  with  increase  in  oxidation  temperature  is 
expected  -  the  oxidation  rate,  and  growth  stress 
i  the  rate  at  which  viscosity  decreases. 


F.  Crystalline  Scale 

Si02  scale  crystallization  kinetics  have  been  measured  for  Hi-Nicalon™-S  SiC  fiber;  crystallization 
begins  after  ~  100  hours  at  1000°C  or  1  hour  at  1300°C,  when  scales  are  roughly  500  nm  and  100  nm 
thick,  respectively.5"6  Unfortunately,  SiC  oxidation  kinetics  for  crystalline  scales  (cristobalite  and 
tridymite)  are  not  well  characterized,  and  there  is  even  less  data  on  creep  rates  and/or  viscosities. 
However,  it  is  clear  that  the  oxidation  kinetics  for  crystalline  scales  are  slower.5'6  Some  kinetic  data 
suggests  that  crystalline  scale  thicknesses  are  ~l/3  those  for  amorphous  scales.55"56  Temperature 
dependence  was  not  reported,  so  to  a  first  approximation  it  is  assumed  that  B0  is  an  order  of  magnitude 
smaller  (B0  =  1.2  x  10"9  m2/s)  and  other  Deal-Grove  parameters  are  unchanged.  Limited  data  for  creep 
rates  for  cristobalite  refractories  suggests  that  at  1550°  -  1650°C  creep  is  negligible  at  very  low  stress 
(0.2  to  0.6  MPa),  but  further  quantification  was  not  possible.36"38 

To  explore  the  possible  effects  of  crystallization  on  growth  stress,  we  assume  that  cristobalite 
viscosity  is  four  orders  of  magnitude  higher  than  that  of  amorphous  silica  at  all  temperatures:  C0  =  3.8  x 
10"9  Pa-s  [24],  Calculations  are  done  for  a  1  pm  thick  scale  at  1100°C,  and  for  cases  where  crystallization 
occurs  when  the  scale  is  100  nm  thick,  500  nm  thick,  and  900  nm  thick  (Fig.  12).  These  thicknesses  set 
the  value  for  X|  in  [1],  Crystallization  rates  are  very  rapid  compared  to  scale  growth  rates.5'6  We  assume 
that  during  crystallization  all  stresses  in  the  amorphous  Si02  are  eliminated  and  are  not  transferred  into 
the  crystalline  product.  Calculations  for  an  Si02  glass  (Fig.  7)  are  included  for  comparison. 

At  1100°C,  Ge  and  gz  are  tensile,  <  200  MPa,  and  in  steady-state  throughout  most  of  the  scale 
(Fig.  12).  For  early  scale  crystallization  (x,  =  100  nm),  Ge  is  1.6  GPa  (tension)  at  the  scale  surface,  and  gz  is 
450  MPa  (tension).  For  x,  =  500  nm,  these  values  are  700  MPa  and  150  MPa,  respectively,  but  a  higher 
tensile  Ge  of  1.5  GPa  is  present  at  x  =  500  nm,  with  a  sharp  transition  to  high  compressive  stress  at 
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Distance  from  SiC-Si02  Interface  (nm) 

Fig.  12.  Hoop  (ere)  and  axial  growth  stress  (ct2) 
levels  after  oxidation  at  1100°C  as  a  function  of 
position  in  a  crystalline  Si02  scale  for  SiC  fibers 
with  6  pm  radius  (b0).  Stress  after  scale 
crystallization  at  x,  =  100,  500,  and  900  nm  is 
calculated,  as  well  as  glass  for  comparison. 


the  point  in  the  scale  at  which  crystallization 
occurred.  For  Xj  =  900  nm  the  scale  has  spent  a 
relatively  short  time  growing  in  the  crystalline 
state.  Surface  stresses  are  compressive,  but  tensile 
a@  of  700  MPa  develops  at  x  =  100  nm. 

As  expected,  much  lower  scale  viscosities 
cause  much  higher  stresses.  Cracks  that  form 
during  oxidation  in  cristobalite  and  tridymite 
scales,  but  not  in  amorphous  scales,5"6  can  clearly 
be  a  consequence  of  the  much  higher  tensile 
surface  stresses  (>  1  GPa)  that  may  develop  in 
crystalline  Si02  scales.  Porosity  present  in  the 
centers  of  crystalline  scales  might  also  be  inferred 
to  be  a  consequence  of  hydrostatic  tension  that 
develops  at  the  point  in  the  scale  at  which 
crystallization  occurred.  Crystalline  Si02  will 
deform  by  dislocation  plasticity  at  high  x; 
mechanisms  implicit  in  the  Eyring  model  for  glass 
no  longer  apply.  There  is  evidence  for  intense 
dislocation  plasticity  in  crystalline  Si02  scale  near 
the  SiC-Si02  interface.5'  The  20  nm  of  Si02 
adjacent  to  SiC  will  have  x  >  1  GPa,  which  is  more 
than  enough  for  domination  of  creep  by 
dislocation-based  mechanisms.57 


compressive  growth  stress  and  externally  ar 
tensile  stress  was  briefly  mentioned.8, 12, 20, 


G.  Effects  of  Stress  on  Oxidation  Rate 

The  retardation  of  silicon  oxidation  by 
plied  stress,  and  the  enhancement  by  externally  applied 
-25, 29-31  50me  data  for  oxidation  enhancement  by  tensile 
stress  also  exists  for  SiC.87  Recent  work  on  SiC  fiber  oxidation  finds  the  oxidation  rates  tend  to  be  higher 
than  those  observed  for  bulk  SiC.5'6  There  was  speculation  that  this  could  be  a  SiC  substrate  geometry 
effect,  related  to  the  development  of  tensile  stress  in  the  fiber  scale. 

The  usual  method  for  handling  stress  effects  on  oxidation  kinetics  is  to  assume  a  pressure  (p) 
effect  on  the  parabolic  rate  constant  B: 


B(p)  =  BePVd/kT  [52] 

where  Vd  is  an  activation  volume,  calculated  to  be  7.5xl0'29  m3  for  silicon  oxidation,12, 20, 25  and  pressure 
is: 

1 

P  =  3  (Ge  +  °z  +  °r')  [53] 

The  pressure  clearly  varies  from  high  negative  values  (hydrostatic  compression)  at  the  SiC-Si02  interface 
to  lower  values  at  the  surface;  these  can  reverse  sign  be  in  hydrostatic  tension  for  thick  scales  at  low 
temperatures  (Fig.  7). 

A  first-order  correction  for  the  relative  effects  of  growth  stress  in  flat-pale  and  cylindrical 
geometries  is  attempted  by  calculating  the  average  self-pressure  (pav)  in  the  Si02  scale  throughout  its 
thickness,  as  a  function  of  total  scale  thickness  (w),  temperature,  and  fiber  radius.  This  is  done  by 
summing  the  pressure  -  scale  volume  product  for  each  scale  increment  and  normalizing  by  the  scale 
volume: 


P(av)  =  Tra  Zj=i  P(bf  -  bjLi)  [54] 

Li  ui 

The  average  pressure  (pav)  was  calculated  for  6  and  3  pm  radius  fibers  (Fig.  13)  as  a  function  of 
total  scale  thickness  from  10  to  3000  nm  and  temperature  from  700°  to  1300°C.  Average  pressures  of 


15 

Approved  for  public  release;  distribution  unlimited. 


1000 


E 

3  300 

Q> 

C 

,u 

iS  100 


30 


3000 


Average  Pressure 
6  pm  Radius 


T(°C) 


3000 


Average  Pressure 
3  um  Radius 


1000 


E 

8 

<v 

c 

.o 

iS 

<v 

s 


300 


100 


30 


Fig.  13.  Average  Si02  scale  self-pressure  (pav)  for  6  and  3  pm  radius  fibers  as  a  function  of  scale 
thickness  (w)  and  temperature 


>2.5  GPa  (compression)  and  >300  MPa  (tension)  were  calculated  for  thin  and  thick  scales,  respectively, 
at  low  temperatures.  Scales  formed  at  high  temperature  do  not  support  significant  self-pressures. 

When  comparing  oxidation  rates  for  bulk  and  fiber  SiC,  the  difference  in  self-pressure  between 
the  two  geometries  is  relevant.  Both  geometries  have  scale  pressure.  Calculations  for  the  flat  plate 
geometry  are  shown  in  figure  14.  The  difference  in  average  pressure  (Apav)  between  a  flat  plate  and  a  6 
pm  radius  fiber  is  shown  in  figure  15.  Fibers  with  thick  scales  formed  at  low  temperatures  have  the 
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Fig.  14.  Average  Si02  scale  self-pressure  (pav)  for 
a  flat  plate  as  a  function  of  scale  thickness  (w) 
and  temperature 
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Fig.  15.  The  difference  in  average  Si02  scale 
self-pressure  (Apav)  between  a  flat-plate 
and  a  6  pm  radius  fiber  as  a  function  of 
scale  thickness  (w)  and  temperature. 
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greatest  Apav,  and  should  therefore  show  the  greatest  enhancement  relative  to  flat-plate  oxidation 
kinetics. 

For  previously  reported  oxidation  data  for  Hi-Nicalon™-S  SiC  fiber,5'6  the  highest  Apav  is  ~500 
MPa  for  ~1  pm  thick  Si02  scales  formed  at  1000°C  (Fig.  15).  By  [52],  the  lower  growth  stresses  (relative 
to  flat  plates)  in  fibers  oxidized  under  these  conditions  would  enhance  B  by  ~8.5x  in  comparison  to  flat 
plates,  increasing  scale  thickness  bv  ~3x.  This  is  at  best  only  a  partial  explanation  for  the  relatively  high 
SiC  oxidation  rates  for  Hi-NicalonT  -S  fiber.  Most  data,  particularly  that  at  700°  -  900°C  and  1200°  - 
1300°C,  has  much  lower  APav  and  the  effects  of  growth  stress  would  be  insignificant. 

It  is  interesting  to  note  that  the  data  of  Gauthier  et  al.  for  stress  dependence  of  the  parabolic 
rate  constant  for  Hi-Nicalon™  SiC  fiber  oxidation  at  800°C  is  consistent  with  an  8x  enhancement  at  ~500 
MPa.32  However,  at  higher  stresses  it  is  no  longer  consistent  with  [52], 32  At  1.276  GPa  there  is  a 
measured  enhancement  of  19x,  but  [52]  predicts  enhancement  by  ~600x.  We  suggest  that  this  is  due  to 
stress  relaxation  is  the  Si02  scale  by  [23],  since  shear  stress  is  well  above  xc. 

IV.  Summary  and  Conclusions 

A  method  to  calculate  the  axial,  hoop,  and  radial  growth  stresses  in  Si02  scales  for  the  2.2x 
volume  expansion  during  SiC  fiber  oxidation  was  developed.  Axial  stresses  are  usually  ignored,  but  for 
structural  fibers  they  have  the  largest  effect  on  strength.  The  method  assumes  oxidation  volume 
expansion  is  three-dimensional  (dilatational)  and  that  the  stresses  resulting  from  constraint  of  that 
expansion  are  relaxed  radially  with  an  Eyring  stress-dependent  Si02  viscosity,  although  other 
appropriate  viscosity  models  can  be  substituted.  The  method  can  be  equally  well  applied  to  fibers  of 
silicon  or  other  materials.  High  compressive  hoop  and  axial  stresses  of  ~25  GPa  are  very  quickly  relaxed 
to  much  lower  values  at  all  temperatures.  Stresses  relax  by  radial  expansion,  which  creates  tensile  hoop 
stress  in  the  outer  scale.  Tensile  hoop  stress  eventually  drives  axial  stress  to  a  tensile  state  by  the 
Poisson  effect.  Tensile  hoop  and  axial  stress  can  reach  values  greater  than  2  GPa  for  long  times  at  700° 
and  800°C.  A  continuum  in  stress  response  from  flat  plate  geometry  to  small  fiber  radius  geometry  is 
demonstrated.  The  accuracy  of  the  growth  stress  calculation  method  is  limited  by  knowledge  of 
accurate  values  for  amorphous  silica  viscosity,  which  can  be  affected  by  incorporation  of  carbon  or  other 
impurities  during  SiC  oxidation. 

Tensile  hoop  and  axial  stresses  reach  steady-state  values  that  can  be  described  by  analytical 
expressions.  At  steady-state,  tensile  stress  relaxation  rates  are  counterbalanced  by  hoop  stress 
generation  tied  to  the  radial  expansion  of  scale  and  the  oxidation  rate.  At  temperatures  over  1200°C  all 
growth  stresses  are  rapidly  relaxed  by  viscous  flow  to  very  low  values.  For  Hi-Nicalon™-S  SiC  fibers,  a 
900°-1100°C  oxidation  temperature  -  300  nm  -  1  pm  scale  thickness  window  may  exist  where  high 
tensile  stresses  at  the  scale  surface  are  likely  to  develop.  These  may  have  negative  effects  on  fiber 
mechanical  properties.  Even  higher  tensile  stresses  may  develop  for  thicker  scales  at  lower 
temperatures,  but  the  times  required  to  reach  those  thicknesses  are  very  long. 

Crystalline  Si02  scales  are  assumed  to  have  much  higher  viscosities  than  amorphous  scales,  and 
should  therefore  develop  much  higher  tensile  hoop  and  axial  stresses  than  those  scales.  Unfortunately, 
lack  of  data  on  cristobalite  and  tridymite  viscosity  or  creep  rates,  as  well  as  diffusion  rates,  prohibits 
accurate  quantitative  prediction.  However,  qualitative  predictions  of  high  tensile  stress  in  the  outer 
scale  and  high  shear  stress  near  the  SiC-Si02  interface  are  consistent  with  microstructure  observations  in 
crystalline  silica  scales. 

Sample  geometry  may  affect  oxidation  rates  through  the  geometric  effect  on  growth  stresses.  A 
Hi-Nicalon™-S  SiC  fiber  and  a  flat  plate  of  the  same  composition  are  calculated  to  have  up  to  a  ~500 
MPa  difference  in  average  Si02  scale  pressure  from  growth  stress  in  thick  scales  formed  at  ~1000°C.  This 
could  increase  scale  thickness  in  the  fiber  by  ~3x  over  that  of  the  flat  plate.  However,  the  effect  is 
smaller  at  higher  and  lower  temperatures,  and  is  currently  only  a  partial  explanation  for  higher  fiber 
oxidation  rates. 
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